Graphene-based resonant-tunneling strucures 
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Resonant electronic transmission through graphene-based double barriers (wells) is studied as a 
function of the incident wave vector, the widths and heights (depths) of the barriers (wells), and 
the separation between them. Resonant features in the transmission result from resonant electron 
states in the wells or hole states in the barriers and strongly influence the ballistic conductance of 
the structures. 
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Recently graphene and graphene-based microstruc- 
tures have been realized experimentally. They exhibit 
unusual properties and have potential technological ap- 
plications especially as an alternative to the current Si- 
based technology. The relativistic-like properties of car- 
riers in grapheneiii^ii^i^, such as an unusual quantum Hall 
effecli^, result from the gapless and approximately linear 
electron spectrum near the Fermi energy at two inequiv- 
alent points of the Brillouin zone. The charge carriers in 
these structures are described as massless, chiral "rela- 
tivistic" fermions, governed by the Dirac equation, and 
their "relativistic" behavior is expected to lead to the ob- 
servation of perfect transmission across potential barri- 
ers, known as Klein tunneling^. In addition, very recently 
electronic confinement was demonstrated in graphene mi- 
crostructures using standard lithography methods^. Re- 
cent studies have considered tunneling through single 
graphene barriers^, wells^, and n-p-n junctions^. One 
important aspect of the electronic transport in quantum 
structures in graphene, that so far has not been investi- 
gated, is resonance effects on the transmission. In this 
letter we investigate them as a function of the barrier 
(well) height (depth) and separation, and contrast those 
for double barriers with those for "non relativistic" elec- 
trons. 

The crystal structure of undoped graphene layers is 
that of a honeycomb lattice of covalent-bond carbon 
atoms. To each carbon atom corresponds a valence elec- 
tron and the structure is composed of two sublattices, 
labelled A and B. The low-energy excitations of the sys- 
tem in the vicinity of the K point and in the presence of 
a potential U are described by the 2D Dirac equation. 
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where the pseudospin matrix a has components given by 
Pauli's matrices, p = (j)x,Py) is the momentum operator. 
The "speed of light" of the system is vp, i.e., the Fermi 
velocity {vp ~ 1 x 10^ m/s). The eigenstates of Eq. 
(1) are two-component spinors ^' = [ipA , ips]'^ , where 
ipA and -tpB are the envelope functions associated with 
the probability amplitudes at the respective sublattice 
sites of the graphene sheet. The term ocmvp creates an 
energy gap and may arise due to an interaction with a 
substrate. 



In the presence of a one-dimensional (ID) confining po- 
tential U = U{x) we attempt solutions of Eq. (1) in the 
form ^A{x,y) = (j)A{x)e'^''^y and ^B{x,y) = i(l)B{x)e^'^^^ 
because of the translational invariance along the y direc- 
tion. The resulting coupled, first-order differential equa- 
tions for 4>a{x) and 4>b{x) can be easily decoupled. They 
read 
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Here ^ = x/L, f7± = e - u ± A, 17 = {VL+n^Y'^, /? = 
kyL, u = UL/hvp, e = EL/hvp, and A = mvpL/h; L 
is the width of the structure, u' the derivative of u with 
respect to ^, and the +(— ) sign refers to C = ^ (C = B). 

For a single quantum well, with a square potential 
of height Uo, Eqs. (2) admit solutions which describe 
electron states confined across the well and propagating 
along it^. For confined states, the spinor components de- 
cay exponentially in the region ^ < —1/2. Then (j)A{x) 
can be written^ as 
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e<-l/2, 

-l/2<C<l/2, (3) 
e > 1/2, 



where a = [{e-uo)"^ - f3^ - A"^]^/^ and = e2_^2_^2^ 
Notice that 0yi(x) depends on ky through (3 = ky/L. A 
similar expression holds for (pB- Then, the transmission 
coefficient is obtained as T = |^3p and ^3 is determined 
by matching 0^ and 0b at ^ = —1/2 and ^ = 1/2. This 
procedure can be repeated for double wells or barriers 
(with K — > — k) to obtain the transmission T through 
them. We avoid giving any expressions for T and instead 
concentrate on the results for it. We obtain them with- 
out the mass term since A is usually very small. 
i) Double barriers. A {ky,kx) contour plot of the loga- 
rithm of the transmission T through a double barrier is 
shown in Fig. 1. Typical values for the barrier height 
Uq — 50 meV, the barrier width L ^ 50 nm, and the 
inter-barrier separation d — 100 nm were used. As seen, 
T depends on the direction of propagation or incident 
angle 9 (tan 6* = ky/kx) and its overall behavior is sim- 
ilar to the earlier single-barrier transmission® with the 
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FIG. 1: Contour plot of the electronic transmission through 
a double barrier for Uo = 50 meV, width L = 50 nm, and 
inter-barrier separation d = f 00 nm. 



exception of well-defined resonances. The overall direc- 
tional dependence of T results from the chiral nature of 
the quasiparticles in graphene and a resonance caused 
by confined hole states in the barrier. Notice the perfect 
transmission T = 1 for normal or near-normal incidence 
(tan 6* « 0), which is a signature of Klein tunneling®. Res- 
onant transmission, which is typical for ordinary double- 
barrier structuresii, is found for values of kx that are 
integer multiples of n/d. 

Another way of presenting the transmission is to plot 
T versus the angle 9 for fixed energy, as shown in Fig. 
2 for £; = 72, 78,85 meV by the red-dashed, black-sohd, 
and blue-dotted curve, respectively. As seen, the trans- 
mission is perfect not only for normal incidence but also 
for particular incidence angles, for which the values of 
the momentum components match the resonance peaks 
shown in Fig. 1. 

It is interesting to explore the transmission as a func- 
tion of the inter-barrier separation d. A (d,ky) contour 
plot is shown in Fig. 3 for Uq = 100 meV, E = 27 meV, 
and L = 50 nm. Notice that for small values of ky, i.e. 
small incidence angles, the transmission is perfect regard- 
less of the magnitude of d. For larger ky the transmis- 
sion oscillates as a function of the inter-barrier separation 
with a period that depends on the momentum. This is 
a new feature, absent from ordinary resonant structures, 
in which T is independent of ky . 

a) Double wells. The result for the transmission of 
electrons, with energies above a single well, was given in 
Ref. 9 and that for a double well is shown in Fig. 4. 
Though some features are qualitatively similar in both 
cases, that for a double well is characterized by a strong 
enhancement of the perfect transmission regions and has 
many resonant features that are absent from that for a 
single well. As compared to the double-barrier case, cf. 
Fig. 1, the resonances are much broader and the trans- 
mission on the average much higher. 

iii) Tunneling current. Selecting the wave vector com- 
ponents kx and ky is possible using quantum point con- 
tacts. However, experimentally one usually measures the 
"average" transmission. That is, one usually measures 
the current J which is proportional to a weighted integral 
of the transmission T{kx,ky). Repeating the argument 



FIG. 2: Transmission through a double barrier as a function 
of the angle of incidence for Uo = 100 meV, L = 50 nm, and 
inter-barrier separation d — 100 nm. 
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FIG. 3: {d,ky) contour plot of the transmission through a 
double barrier for Uo = 100 meV, E = 27 meV, and L = 50 
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FIG. 4: As in Fig. 1 but for a double well with 
energy E > Uo, Uo — 50 meV, and L — 200 nm. 




FIG. 5: Conductance through a single and double barrier vs 
energy E for A = 0. Inset; conductance through the same 
structures for mvp — 270 meV; the red (green) curve is for a 
single (double) barrier. 
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FIG. 6: Conductance through a single and double well vs 
Fermi energy Ef- The parameters are given in the text. 

of Ref. 11 for 3D electrons with a parabolic energy spec- 
trum, we find that for the linear spectrum E = Hvph the 
current J, due to a voltage drop eV along the x direction, 
is given by (A = 2evF/h?) 

J = -\J T{E,9)[f{E)- f{E + eV)]EdE cos 9d9, (4) 

where f{E) is the Fermi function. For low temperatures 
we can approximate [f{E) — f{E+eV)] by —eV5{E—EF) 
and extract the low-bias conductance from Eq. (4). For 
a 2D parabolic spectrum A is replaced by 2e^/2m/h? and 
EdE hy E'^/^dEiYiY.a^. (5). The main difference between 
the two cases is that in the latter the transmission T{E, 6) 
depends only on the magnitude of the a;-component of the 
momentum, whereas in the former it depends on both 
components. 

In Fig. 5 we plot the conductance vs the Fermi en- 
ergy for a single-barrier (dashed curve) and for a double- 
barrier structure for symmetric (solid curve) and asym- 
metric (dotted curve) barriers with width L = 50 nm. 
The height of the first barrier is C/q = 100 meV and that 
of the second Uq — 100, 50, and meV, for the solid, 
dashed, and dotted curve, respectively. The inter-barrier 
separation is d = 50 nm. For energies below the max- 
imum barrier height, the transmission is dominated by 
the Klein tunneling effect, in which the incident elec- 
trons are resonantly transmitted via the confined hole 
states in the barriers, resulting in peaks in the conduc- 



tance. An additional pronounced structure is also found 
for double barriers due to an additional resonance effect 
caused by the quasi-confined electron states in the well 
region. This resonant structure is most explicit for sym- 
metric barriers. One striking consequence of that is the 
fact that the low-bias conductance can be significant low- 
ered as the height of one of the barriers is decreased. As 
the Fermi energy increases above the maximum barrier 
height (here 100 meV ) the resonant features practically 
disappear. In contrast with the A ^ results shown in 
the inset, which are a good qualitative approximation of 
the parabolic spectrum case, the results for graphene, in 
the limits Ep — > and Ep — > c», show T — 1 for near- 
normal incidence, see, e.g.. Fig. 1. Using Eq. (4) we 
find immediately G oc Ep which agrees with our numer- 
ical results of Fig. 5. The present results are valid only 
in the ballistic regime, but in the presence of disorder 
the results for the conductance near the Dirac point is 
expected to be strongly modified. 

In Fig. 6 we plot the conductance vs the Fermi energy 
for a single well (dashed curve) and for a double well, 
symmetric (solid curve) and asymmetric (dotted curve). 
The well width is 50 nm. The depth of the first well is 
Uq = 100 meV and that of the second C/q = 100, 50, and 
meV, for the solid, dashed, and dotted curve, respec- 
tively. The inter-well separation is c? = 50 nm. In this 
case, only weak resonant features are found that are de- 
termined completely by the first well. The presence of the 
second well almost does not influence the conductance. 

In summary, we reported unusual resonant-tunneling 
features in graphene microstructures and a strong direc- 
tional character of the transmission through them. The 
latter is perfect for a range of incident angles. No such 
behavior is found in the case of a parabolic spectrum. 
The conductance displays several resonant features that 
can be tuned by the barrier heights as well as by the 
separation between the barriers. 
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